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Abstract 

The effect of Majorana neutrinos in the MeV mass range on the double beta decay of various 
isotopes is studied on pure phenomenological arguments. By using only experimental half life data, 
limits on the mixing parameter U^ h of the order 10~ 7 can be derived. Also the possible achievements 
of upcoming experiments and some consequences are outlined. 

One of the major challenges of modern physics is the still open question of a nonzero neutrino mass. 
A massive neutrino could have important consequences for astrophysics and cosmology, i.e. massive 
neutrinos are good dark matter candidates and massive neutrinos are the preferred solution of the solar 
neutrino problem via the MSW-effect. From the particle physics point of view massive neutrinos would 
require a modification of the successful standard model and would open one of the most promising ways 
for testing models beyond the standard model (GUT-theories). For the physical potential of massive 
neutrinos see (lj. 

Direct measurements of neutrino masses result at present in limits for the three neutrinos of 

m„ e < 4.35 eV (out of tritium beta decays [2] ) , (1) 
m v < 160 keV (out of pion-decays || ) (2) 
m„ T < 24 MeV (out of tau-decays |Q) (3) 

As can be seen MeV r-neutrinos are not ruled out at present. In recent times there is a growing interest 
in models with MeV r-neutrinos ||. To be cosmological acceptable such neutrinos must be unstable 
because otherwise they will overdose the universe. On the other hand bounds on the number of neutrino 
flavours coming out of big bang nucleosynthesis have been relaxed recently, now allowing a value between 
about 2.2 @ and 3.9 §. This opens space for a MeV-mass of v r . 

In this paper the effects of a MeV Majorana-neutrino in double beta decay are investigated. For a 
discussion of heavy sterile neutrinos in double beta decay see The analysis follows partly that of 
H [To). Neutrinoless double beta decay {Quj3j3 decay ) of a nucleus (A,Z) 

(A,Z) (A,Z + 2) + 2e- (4) 

violates lepton number by 2 units and is only possible if neutrinos are Majorana particles (see Q). There 
are about 35 possible double beta decay emitters, about 10 of them have experimental obtained half life 
limits of larger than 10 20 y. At present the best limit results from the Heidelberg-Moscow collaboration 
studying the decay of 76 Ge [jllj 

T°/ 2 > 7.4 • 10 24 a -> < m„ c >< 0.56 eV (90% CL) (5) 
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The measured quantity (neglecting right handed weak currents), called effective Majorana mass < to„ c >, 
is given in the case of light neutrinos (m, < 1 MeV) by 

jv 

<m Ve >=\Y,U! l m l | (6) 

i=l 

where rrii characterizes the N mass eigenstates and J7 e , the mixing matrix elements. Things change in 
case of heavy neutrinos (m v > 1 MeV). By calculating the nuclear matrix elements for double beta decay 
involving MeV neutrinos the neutrino mass in the neutrino propagator can no longer be neglected with 
respect to the neutrino momentum. For a detailed discussion on the matrix element calculations see jl2| . 
This results in a change of the radial shape of the used neutrino potential H(r) from 

, . 1 . . . exp(— m^r) . . 

H (r) oc — light neutrinos — > H (r) oc heavy neutrinos (7) 

r r 

This changes can be accomodated for by introducing an additional factor F(rrih, A) in eq. (|^), which 
depends on the mass of the heavy neutrino and on the atomic number A of the nucleus. Eq. (j(|) is 
modified to 

N M 

<m Ve >=\ U li m i+ F{m h ,A)U 2 eh m h \ (8) 

i—l,light h—1, heavy 

Assuming one heavy neutrino with nih — rn 2 the function F(m,h, A) is given by 

F(m 2 ,A) =< 1/r > _1 < exp(-m 2 r)/r > (9) 

r corresponds to the distance of the two nucleons in the nucleus undergoing 0v(3j3 decay . The average 
is with respect to the two nucleon correlation function appropriate for the nucleus. Using a correlation 
function containing a hard core repulsion characterized by a hard core radius r c of 0.5 fm between the 
two nucleons 

= 4/3^)3 - rg] g ( f -^(2fl-r) (10) 

it follows 

F(m 2 ,A) = 0.5- — r[(l + m 2 r c )exp(-m 2 r c ) - (1 + 2m 2 R)exp{-2m 2 R)} (11) 

(m 2 Ry 

For all nuclei of interest (A > 48) the nuclear radius R is much larger than the hard core radius (R>r c ). 
Therefore ^(7712, ^4) varies with the nuclear radius as R~ 2 or using the relation R ~ 1.2A 1 / 3 fm it results 
in a dependence of A~ 2 / 3 . 

Consider now the simple case of an electron coupled via the standard weak charged current to two massive 
Majorana neutrinos xi,2 under the assumption of CP-conservation: 

v e = xicosO + X2sin9 (12) 

where the fields Xi,2 satisfy the Majorana condition (C is the charge conjugation matrix): 

Xi,2 = 171,2^x1,2 Vi = ±:L ( 13 ^ 
The corresponding < m„ e > is then given as 

< m Ue > —\ m\cos 2 Q + rjirj 2 m 2 sin 2 9 | (14) 
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in the case of only light neutrinos or for one light and a MeV-neutrino as 

< m Ve > =| micos 2 9 + F(rri2, A)i]ir)2'm2sin 2 9 | (15) 

In case of small mixing angles 9 the /3-decay experiments measure mi. For double beta decay if 9 ^ 
and rj\rj2 = -1 something like destructive interference can occur in < m„ e > and the value measured can 
be smaller than the one from the tritium experiments. 




Figure 1: left: Limit on the mixing angle sin 2 9 e i as a function of the heavy neutrino mass m,2 in the 
region of 1-100 MeV for five combinations of isotopes. At present the best limit is given by the Ge-Te 
pair. The combinations of the Ca-Nd and Ca-Te result are nearly identical, right: Upper limit of the 
Majorana mass of a light neutrino mi in the case of interference with one heavy neutrino. The Ge-Te 
pair results in the case 0/7712 ~ WOMeV in a limit for mi of about 6 eV. 

Using eq.(|il"|) and ( |l5| ) it can be seen that < m Ve > gets an A-dependence because of F(m,2, A) making it 
worthwile to look into experimental results of different double beta decay emitters. In Tab.l a comparison 
of some double beta decay emitters as well as present limits on the half life T®" 2 and the effective Majorana 
neutrinos < m v >are shown. 



Isotope 


4 «Ca 


/6 Ge 


m Mo 




m Xe 




Tty 2 present 


9.5 E21 


7.4 E24 


4.4 E22 


7.7E24 


4E23 


2.1E21 


< m Ue > 


12.8 


0.56 


5.4 


1.0 


2.4 


4.0 


T{% future 


1 E23 


1.5 E25 


1E25 


7.7 E24 


3.6 E25 


1E23 



Tab.l: Some selected isotopes with reasonable half life limits for double beta decay . Shown are the present 
limits for the half life T^ 2 , the resulting < m„ e >- limit using the matrix elements from JTjj/ , as well as 

some proposed limits of ongoing, upcoming or planned experiments on T®Y 2 . 
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Clearly there are two strategies to follow for looking at the A-dependence: The largest effect is expected 
if the two isotopes have the largest possible spread within A f.e. 48 Ca and 150 Nd or even 23S U. On 
the other hand much better experimental limits exist for the 0v/3(3 decay half life for isotopes like 76 Ge 
and 128 Te . Fig. la gives an idea of the variation in the mixing angle as a function of the heavy neutrino 
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Figure 2: Comparison of the Ge-Te double beta decay limit (solid line) on the mixing angle with other 
experiments. Shown are the results from [15] (open circles), [16] (asterix), [17] (stars) and [18] (filled 
circles). It can be seen that the double beta decay results give the best limit below about 35 MeV. The 
region above the curves are excluded. 



76 Ge (best half life limit) and 
Nd (largest A with reasonable 



mass for five different pairs of isotopes. On the "small" A side it is 
is Ca (lowest A) on the "higher" side 128 Te (best half life limit) and 150 
half life limit). It can be seen that the much better experimental half life limits of 76 Ge and 128 Te 
overcompensate the smaller separation in A. Therefore this pair will be used for a comparison with other 
experiments. Alternatively we can restrict the limit on the light neutrino mass mi assuming a nearly 



perfect cancellation of < m Vp > in 76 Ge by using a relation between two isotopes like [[i~4| 



mi < 



F(m 2 ,76) 



F(m 2 ,76)-F(m 2 ,A) 



-(< m„(76) > + < m v (A) >) 



(16) 



where A is one of the other isotopes. The value of the observable light neutrino mass as a function of 
the heavy one is shown in fig. lb. Using the simple mixing scheme of eq.(|l2|) fig. 2 shows the limits on the 
mixing parameter | U e i \ 2 — sin 2 9 e i because of the A-dependence in double beta decay in comparison with 
other experiments. It can be seen that in the region below about 35 MeV this limit is the most stringent 
one. This can be improved by the ongoing, upcoming or planned future experiments. Using the proposed 
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Figure 3: Same presentations as fig.l for ongoing, upcoming or planned experiments and their proposed 
half lifes limits. The best mixing angle limit (left) is now given by the Ge-Xe pair, which also gives an 
upper bound on m\ of 2 eV assuming mi re lOOMeV. 



half live limits of the experiments shown in Tab.l limits on the mass eigenstate and the mixing angle can 
be obtained as shown in fig. 3. In this case the 76 Ge and 136 Xe data would dominate the bound, which 
is by a factor of about 3 lower than the present limits. So far all the considerations are based on very 
general arguments. 

Now look at the standard model and assume that there are no extra neutrinos from yet unknown physics. 
The Z°-resonance width results in 2.983 ± 0.025 Jl9| neutrinos flavours lighter than 45 GeV. Using the 
upper bounds on the different neutrino masses as given in eq. (1-3) it is obvious that m<i can only be 
associated with v T . Therefore the derived values for | t/ ei | 2 correspond to | U eT 2 for Am 2 re 10 12 eV 2 . It 
should be mentioned that the bounds on | U eT | 2 are weakened by some orders of magnitude by including 

and assuming that there is a cancellation of the and v r contributions to < m Ui! >. 
A direct consequence of the derived limit can be applied to the decay of v T — > v e + e + e~ because of 
the lifetime dependence of r oc| U e i | 2 (neglecting other decay channels, see ||). Different limits on the 
lifetime of v r as a function of its mass including the one obtained via double beta decay can be seen in 
fig. |i[ By using only rather model independent limits most of the allowed parameter space can be ruled 
out. 

As a second example consider the possible evidence for neutrino oscillations as seen with the LSND- 
detector ^J. This would require a massive neutrino of re 1 eV. Assume that the mass of the electron 
neutrino is somewhere in the region of 1-4 eV, which is still allowed by beta decay experiments, a heavy 
neutrino between 1-24 MeV, and that a perfect cancellation in 76 Ge Qv(3j3 decay occurs. This would 
imply values of < m Ve >in other isotopes like 100 Mo, 128 Te and 150 Nd in the range 

2.7 • 10~ 3 ey <<m Vs >< 0.25eV ( 100 Mo) (17) 
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Figure 4: Lifetime of v T — > ^ e + e + e~ against the tan-neutrino mass ra Vr . Shown is the exclusion due to 
double beta decay (lower curve) The vertical line at 24 MeV corresponds to the ALEPH-limit, the upper 
curve is excluded by cosmological arguments (Overdensity of the universe). These two bounds are valid 
for any decay mode. Also shown is the bound on the v T — > v e M +7(e+e~) decay-mode given by SN 1987a 
[20]. ' 

5.5 • lCT 3 eU <<m Ue >< 0.49eU ( 128 Te ) (18) 
7.4 • l(T 3 eU <<m„ e >< 0.63eU ( 150 Nd) (19) 

The upper bounds are in the region of upcoming or planned experiments and such a scenario can be 
tested within the near future. 

Conclusions 

The effects of MeV Majorana neutrinos in double beta decay are investigated by comparing half life limits 
of different isotopes. Using only experimental bounds the most stringent limits on the mixing matrix 
element U e i in the region from 1-35 MeV are obtained. Using the bound on U e i a large part in the 
t„ t — m„ T plane for the decay mode v T — ► v e e + e~ can be excluded. Whether there is a Majorana v e 
somewhere in the eV region and a MeV v T can be tested or detected within the near future. 
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